Absence of gross inhomogeneity in HTS cuprates 
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Recent scanning tunnelling microscopy (STM) studies on Bi2Sr2CaCu20s+5 (Bi-2212) revealed 
the presence of severe inhomogeneity with length scale Lo ~ £o, the coherence length. Other 
studies have been interpreted in terms of mesoscale or nanoscale phase segregation. Here we analyze 
heat capacity and NMR data on Bi-2212 and (Y,Ca)Ba2Cu307_<5 and find no evidence for phase 
segregation or gross inhomogeneity. For Bi-2212 the coherence scale Lo increases with doping from 
5 to 17£o and the hole density inhomogeneity decreases from 0.028 to 0.005 holes/Cu. We conclude 
that STM measurements considerably overstate the inhomogeneity in bulk Bi-2212. 
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The possible importance of inhomogeneity in the 
physics of high temperature superconducting (HTS) 
cuprates has been stressed since the earliest investiga- 
tions of these materials. Of necessity they involve chemi- 
cal disorder as substitutional dopants in the parent com- 
pounds are usually required to introduce carriers. Fre- 
quently observed forms of inhomogeneity include substi- 
tutional, interstitial, reconstructive, twinning and incom- 
mensuration [1—4]. Despite this potentially high degree 
of structural and compositional disorder (which is quite 
variable from one compound to another) the cuprates 
display universal phase behaviour [5-7]. Their physi- 
cal properties vary systematically with doping in very 
much the same way from one compound to another pro- 
vided simply that the substitutional disorder lies out- 
side of the active Cu02 planes. This, of course, reflects 
the strongly two-dimensional nature of the electron dy- 
namics in the CuC>2 planes. More fundamentally how- 
ever, inhomogeneity of intrinsic electronic origin has been 
discussed widely, most notably various nanoscale stripe 
phases in which the background spins separate spatially 
from doped charges to form quasi-one-dimensional struc- 
tures [8] . There is thus much scope in HTS materials for 
significant extrinsic and intrinsic spatial inhomogeneity. 

Recent high-resolution scanning tunnelling microscope 
(STM) studies [9,10] at 4.2K have provided dramatic ev- 
idence of electronic inhomogeneity on the surfaces of Bi- 
2212 in the superconducting (SC) state. These reveal a 
patchwork of contrasting regions with length scale w 14A, 
comparable with the coherence length £o, some having 
well-defined coherence peaks and others having strongly- 
reduced coherence peaks shifted to higher energies similar 
to the features of the pseudogap above T c [11]. The ~ 
50% spread in the magnitude of the SC energy gap, A , 
was shown to correlate with a similar spread in the local 
density of states. These variations were attributed to a 
distribution of hole concentration, p, arising from oxygen 
disorder and a FWHM spread of Ap w 0.08 holes/Cu was 
inferred by comparing the distribution in SC gap with 
the p-dependence of A deduced from ARPES. It was 



suggested that such inhomogeneity may be quite general 
and may reflect " an intimate relationship with supercon- 
ductivity" in the cuprates [9]. However, it cannot be pre- 
sumed that such inhomogeneity is generic nor that it even 
reflects the nature of the bulk electronic state. STM is 
strictly a surface probe susceptible to just the outermost 
CuC>2 layer. Little is known about surface reconstruction 
in this material, the location of additional oxygen atoms 
in the cleaved BiO layer or the role of Bi atoms incorpo- 
rated into the underlying SrO layer. We argue here from 
thermodynamic measurements and NMR data that the 
bulk electronic state is much more homogeneous than is 
suggested by the STM results which, by implication, are 
an artifact of a perturbed surface. 

A quantitative measure of the degree of doping in- 
homogeneity in the HTS cuprates can be inferred from 
NMR linewidths. NMR is a bulk probe which is sensitive 
to the local electronic state near the probe nucleus. We 
show in Fig. 1 the 89 Y Knight shift, 89 K s , at 300K for 
Y .8Cao.2Ba2Cu3C , 7_ ( s (Y,Ca-123) at six different values 
of S and p spanning the SC phase diagram [12]. The p 
values were obtained from thcrmopower measurements 
and T c values are also plotted. As is well known, the 
Knight shift varies systematically with doping and the 
figure reveals a slope of 580 ppm/hole. Doping inho- 
mogeneity will therefore induce a minimum line width 
of 580 xAp ppm. Other contributions, such as local 
moments [13] and the effect of disorder on the chemical 
shift will also add to the linewidth. We may thus use the 
observed linewidth to determine the maximum doping 
inhomogeneity (if all the linewidth came only from Ap). 

The vertical bands in the figure show the FWHM in the 
hole concentration determined in this way [14]. We have 
used the magic-angle-spinning data of Balakrishnan et al. 
[15] for YBaaCusOr^ (Y-123) and of Williams et al. [16] 
for YBa2Cri408 (Y-124). In the last case the linewidth 
is extremely narrow, about 100 Hz, giving Ap < 0.0069. 
It is clear that this class of HTS is remarkably homoge- 
neous yet they exhibit the same generic thermodynamic 
and ground-state behavior as Bi-2212 [17]. 
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FIG. 1. The doping dependence of the 89 Y Knight shift 
at 300K for Yo.8Cao.2Ba2Cu 3 7 -«. From this the NMR 
line-width is used to estimate a maximum FWHM, Ap, in 
doping state shown for Y-123 and Y-124 by the vertical bands. 

Heat capacity studies on HTS cuprates provide a num- 
ber of potential checks against electronic inhomogeneity. 
These include (i) the value of the specific heat coefficient, 
7 = Cp/T at T—0, (ii) its value, -f n , at high tempera- 
ture in the normal state, (iii) the magnitude and dop- 
ing dependence of the jump, A7, in 7 at T c and (iv) 
the transition width. For example, mesoscale segrega- 
tion into SC and normal metallic phases will result in 
normal states in the gap manifested by a non-zero value 
of 7 at T=0. Progressive segregation with overdoping, 
for example, will see 7(0) rising steadily, filling out the 
d-wave gap. On the other hand, mesoscale segregation 
into insulating and normal metallic phases will result in 
a progressive reduction in the high-temperature value of 
7„ as states are removed. Here we focus on the thermo- 
dynamic data for Bi-2212 because this is the system in 
which the STM inhomogeneity was observed. To this we 
add data for Y,Ca-123, a very different system, to show 
that our observations are generic to the cuprates. 

Ceramic samples (« 2 gm) of Bi2.iSri.gCaCu208+c5 
[18] and Y .8Cao.2Ba2Cu 3 07_5 [19] were synthesized us- 
ing standard solid-state reaction and were oxygen loaded 
by slow cooling and annealing for many days at 325°C in 
flowing oxygen. The electronic heat capacity was mea- 
sured using a high-precision differential technique de- 
scribed elsewhere [20] . The samples were then alternately 
annealed and quenched into liquid nitrogen and the heat 
capacity remeasured so as to track the thermodynamic 
state as a function of progressively depleted oxygen con- 
tent. Changes in oxygen stoichiometry were determined 
from the mass changes and these concurred rather pre- 
cisely with changes in the phonon specific heat. Values of 
p were determined from T c values using the approximate 
parabolic phase curve, T c (p), given by [5] 

T C ( P ) - T c , max [ 1 - 82.6 (p - 0.16) 2 ]. (1) 

Similar studies have correlated such values of p with the 
thermopower for both Bi-2212 [21] and (Y,Ca)-123 [22]. 

Fig. 2 shows a series of 7(T) curves for each sam- 
ple at closely spaced p- values from deeply underdoped to 



deeply overdoped. Dashed curves denote optimal doping 
where T c maximises while the thick curves denote crit- 
ical doping (p crit =0.19) where the pseudogap abruptly 
closes. Apart from small impurity-induced upturns for 
T < 10 K 7(0) remains close to zero across the entire 
doping range 0.084<p<0.23 in Y-123 and 0.13<p<0.21 
in Bi-2212, while the values of 7(T) at high temperature 
are constant, independent of doping and characteristic of 
a uniform normal metal. These observations rule out the 
two types of mesoscale segregation discussed above. We 
note that the 17 Knight shift in overdoped Y,Ca-123 
also approaches zero at T—0 [12] in agreement with the 
thermodynamic data and independently confirming the 
absence of unpaired quasiparticles. Significant gap filling 
in the single layer cuprates Tl-2201 [23] and La-214 [24] 
only occurs in the most heavily overdoped region p>0.24. 
This may reflect the increased role of pair breaking rel- 
ative to a rapidly decreasing superconducting gap, and 
does not necessarily indicate segregation. 
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FIG. 2. The specific heat coefficient, 7(T), for a series 
of doping states for Bi-2212 (0.13 < p < 0.21) [17] and 
(Y,Ca)-123 (0.084 < p < 0.23) [18]. Arrows indicate increas- 
ing doping, bold curves denote critical doping, p = 0.19 and 
dashed curves optimal doping, p — 0.16. 

Fig. 2 also reveals rather sudden changes in j(T,p) 
with doping. The jump, A7, at the transition changes 
only slowly across the deeply overdoped region but, start- 
ing at p cr it=0.19, the magnitude of A7 falls very rapidly 
with further reduction in doping. Adjacent curves are 
separated by about Ap « 0.008 and the observed sud- 
den collapse of the specific heat anomaly suggests that 
any spread in doping state arising from nanoscale inho- 
mogeneity must be less than this value. Consistent with 
this inference the transition widths remain quite narrow 
and we now focus on this issue in more detail. 

A finite transition width At m may arise from a spread, 
Ap, in doping over mesoscopic regions (due for example 
to a spread, Aa;, in local donor density, x). Here, 

At m = AT C /T C = T~ x I dT c /dp I xAp, (2) 
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which falls to a minimum at optimum doping p = p opt 
= 0.16. On the other hand a finite transition width 
Atf s independent of | dT c /dp | will result from inhomo- 
geneity on a much shorter length scale, £ the coherence 
length, due to finite — size effects [25-27]. In the ther- 
modynamic limit the coherence length, £(T), diverges as 
t = T/T c - 1 -> according to the power law £(T) = 
&/ \t T, where ^=0.67 (« 2/3) for 3D-XY fluctuations. 
Where the system is inhomogeneous this divergence is 
cut off by a length scale i which is presumed by several 
authors [25-27] to be the spatial length scale of the in- 
homogeneity. As a consequence the divergence of critical 
fluctuations is suppressed within 



At fs = (£o/£o) 3/2 . 



(3) 



It is immediately apparent that, if L w £ as inferred 
from STM [9,10], then the transition broadening must 
be comparable to the magnitude of T c , which is clearly 
not the case. In fact, typical transition widths of « 2-3K 
require that £o/£o > 10. For any doping state we can 
estimate the transition width from the finite cut-off of 
the fluctuation- induced heat capacity as | i |— > and we 
will see that it is indeed small. 
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FIG. 3. (a) The specific heat anomaly near T c for Bi-2212 
and YBa2Cu30e.92, both with p ~ 0.18. Panels (b)and (c) 
show AC versus log 10 \t | (solid curve) and versus log 10 (t*) 
(dashed curve), where t* = (t 2 + At 2 ) 1 '" 2 . Single-head arrows 
show the transition width, At, and double-head arrows show 
the mean-field jump in specific heat, AC MF . 

Fig. 3(a) shows the specific heat near T c of Bi-2212 
and pure YBa2Cu3 6.92 7 both at a similar doping state 
p « 0.18. In the 3D-XY model the logarithmic divergence 
of the specific heat near T c is given by AC P = A~ In | t \ 
+ ACp 4F for t <0 and A+ In | t | for t >0. Here AC p MF 
is the mean field step, A~ w A + = 4/c B /(9"7r 2 £l(0)) and ft 
is the coherence volume [28]. A semilog plot of AC P ver- 



sus In |t | gives two parallel lines offset by the mean- field 
step, AC p MF . In practice this plot (shown in Fig. 3(b) 
for Bi-2212 and in Fig. 3(c) for Y-123) exhibits negative 
curvature for sufficiently small \t | due to the finite-size 
cut-off in critical fluctuations or other sources of broad- 
ening. The effect of a spread in T c may be modelled by 
replacing t by t* = (t 2 + At 2 ) 1 / 2 in the above expressions 
for AC P [26]. This is illustrated in Figs. 3(b) and (c) by 
selecting values of the half-width At that just avert nega- 
tive curvature. Since Atf s and At m < At the measured 
broadening At yields a lower limit to L /£ and an upper 
limit to Ap (eqs. 2,3). From the figures one finds At = 
0.018 and AC* fF = 34 mJ/g.at.K for Bi-2212; and At = 
0.0028 and AC p MF = 266 mJ/g.at.K for YBa 2 Cu 3 6 .9 2 . 
These values for At translate to Lo > 16 £o and 50 £0, re- 
spectively, values which are very much larger than those 
implied by the STM results. 

hole concentration, p 




FIG. 4. The half-width of the transition At and lower lim- 
its for Lo/£o derived from eq. (3) for (a) Bi-2212 and (b) 
(Y,Ca)-123. Values of At m calculated using eq. (2) are shown 
in (a) for Ap = 0.005 and in (b) for Ax = 0.01. 

The same analysis may be carried out for each inves- 
tigated doping level in Y .8Cao.2Ba2Cu307_5 and in Bi- 
2212. Fig. 4 shows the measured broadening At and 
the corresponding lower limit to L / Co plotted as a func- 
tion of p in the case of Bi-2212 and of oxygen content, 
x = 1 — S, for (Y,Ca)-123. In the former case lower-limit 
values of Lo/Co range from 17 at critical doping (p cr it 
= 0.19) down to 5 in the most heavily underdoped sam- 
ples. For (Y,Ca)-123, though less than for pure Y-123, 
values of L /£ are still as high as 35. Even in the most 
extreme cases the length scale is almost an order of mag- 
nitude longer than that inferred from the STM studies. 
We also show in Fig. 4 At m obtained from eq. (2) using 
Ap = 0.005 for Bi-2212 and Ax = 0.01 for (Y,Ca)-123. 
In the latter case At™ follows At rather well across the 
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phase diagram indicating a spread in oxygen content of 
Ax w ±0.01 for this compound. In Bi-2212 the transi- 
tion width At does not correlate with At m , confirming 
the existence of short-length-scale inhomogeneity. Nev- 
ertheless we can conclude that the spread in p does not 
exceed « 0.005 at p= 0.21 and w 0.028 at p=0.13, signif- 
icantly less than the value inferred from the STM data. 

We plot in Fig. 5 the magnitude of the deduced mean- 
field step, ACp IF , for Bi-2212. This remains rather con- 
stant on the overdoped side but falls extremely sharply 
at Pcrit = 0.19 as the pseudogap opens and, remarkably, 
is essentially zero by optimal doping. A mere change of 
Ap = 0.03 results in the complete collapse to zero of the 
mean-field step. This abrupt change again underscores 
the degree of bulk electronic homogeneity that must be 
present, even in the case of Bi-2212. With AC p MJ? « 
for p < 0.16 the specific heat anomaly becomes symmet- 
ric about T c because its weight arises purely from critical 
fluctuations. We note that there is no significance in the 
fact that the mean- field step reaches zero at p op t- This 
is not the case in Y-123 and we expect that the collapse 
in AC^ F would be even more abrupt in a system with 
weaker interlayer coupling than Bi-2212. The coherence 
volume, O(0), is also shown in the figure. The values are 
consistent with £ b w 17 A and £ c w 0.5 A. 
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FIG. 5. The doping dependence of the mean field jump, 
AC MF , in specific heat for Bi2.iSri.9CaCu20s+5 determined 
from the fluctuation analysis as in Fig. 2(b) for each anneal- 
ing state. The coherence volume 51(0) from the amplitude 
A+ = 4fc s /(97r 2 fi(0)) is also shown. 

We conclude that specific heat and NMR measure- 
ments indicate levels of inhomogeneity much lower than 
those inferred from STM experiments. In the dop- 
ing range considered there are no phase-separated "nor- 
mal" or "insulating" regions coexisting with supercon- 
ductivity. Y-123 and Bi-2212 reveal very narrow ther- 
modynamic features, such as the jump in 7 = C p /T, 
which change abruptly with doping indicating a spread 
in local hole concentration no more than Ap « ±0.01. 
Such features, found also in Yi_ a; Ca a ;Ba2Cu307_ ( s, 
Tlo^Pbo^S^Cai-aYrcC^Oy and in I^-xSr^CuO^ are 
generic to the cuprates and suggest (at least for p >0.125) 
a remarkable degree of bulk electronic homogeneity in 



spite of considerable disorder in the non-Cu02 layers. 
Analysis of critical fluctuations near T c indicates a spread 
of local donor densities no greater than Ax w ±0.01 for 
Yi_ x Ca x Ba 2 Cu 3 07-,5 and w ±0.02 for Bi-2212, while 
the latter compound exhibits a sudden and complete col- 
lapse of the mean-field step between critical and optimal 
doping. We conclude that the static inhomogeneity ob- 
served in STM on Bi-2212 single crystals is unrelated to 
bulk superconducting HTS properties and probably re- 
flects the presence of disorder at the surface. 
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